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ABSTRACT 

Cool stars like our Sun are surrounded by a million degree hot outer atmosphere, the corona. Since 
more than 60 years the physical nature of the processes heating the corona to temperatures well in 
excess of those on the stellar surface remain puzzling. Recent progress in observational techniques and 
numerical modeling now opens a new window to approach this problem. We present the first coronal 
emission line spectra synthesized from three-dimensional numerical models describing the evolution 
of the dynamics and energetics as well as of the magnetic field in the corona. In these models the 
corona is heated through motions on the stellar surface that lead to a braiding of magnetic field lines 
inducing currents which are finally dissipated. These forward models enable us to synthesize observed 
properties like (average) emission line Doppler shifts or emission measures in the outer atmosphere, 
which until now have not been understood theoretically, even though many suggestions have been 
made in the past. As our model passes these observational tests, we conclude that the flux braiding 
mechanism is a prime candidate for being the dominant heating process of the magnetically closed 
corona of the Sun and solar- like stars. 

Subject headings: Sun: corona — stars: coronae — Sun: UV radiation — MHD 



1. INTRODUCTION 

Shortly after it was realized in the 1930 ies that the 
corona is hot ijGrotrianI [19391 lEdlenl 1X943) . first pro- 
posals for the heating mechanism were made based on 
upward pro pagating sound wav es generated on the stel- 
lar surface l|Schwarzschildl[T948|l because convective mo- 
tions contain far more energy than what is needed to 
heat the corona. Later it became clear that the heat- 
ing mechanism has to be related to the magnetic field 
dominating in the corona. The coronal magnetic field 
is rooted in the photosphere, however, where it is dom- 
inated by the kinetic energy of the convection. Thus 
at the surface magnetic field-lines are pushed around, 
and the (stochastic) footpoint motions result in a braid- 
ing of the magnetic field lines in the corona. This 
implies that magnetic field gradients are built up, in- 
ducing currents which are finally dissipa t ed and thus 
heat the co rona as suggested by IParked lll972L I1994D 
— see also Sturrock & Uchida fl981"): 'Parker fl983'); 
Ivan Ba llegooiicn (1986); Hcvvacrts & Priest (.1992). 
Thus the magnetic field acts as the agent to channel 
energy from the cool surface into the hot corona. Sev- 
eral studies have been carried out to investigate the role 
of footpoint motions for coronal heating. For example 
IPricst & Schriivcr (1999) studied the influence of moving 
magnetic sources of opposite polarity in the photosphere 
on reconnection in the above corona and Priest et al.l 
((2002) constructed a model view for coronal heating 
based on flux-tube tectonics related to the ever changing 
structure of the magnetic carpet (Schriivcr et al. 199^. 
Many other mechanisms have been discussed to heat the 



corona, e.g. based on waves, but there are too many of 
them to be reviewed here ( see the S OHO 15 proceedings 
for an up to date overview. iDanesvl 12^0 041 . 

Even though the proposal of field line braiding has been 
widely discussed, it was not until recently that numeri- 
cal modeling made it possible to test whether this mecha- 
nism actually works. Studies had been conducted for the 
flux braiding on small scales, investigati ng the amount 
and character of the energy depo sition l|Hendrix et alJ 
119961 TCalsgaard & Nordlundl [T99^ . Finally a smafl ac- 
tive region could be described in a numerical experiment, 
resulting in a realistic looki ng corona consisting of a num- 
ber of loop-like structures l|Gudiksen fc Nordhmdll2?)(l^ 
I20n4blan . 

In this work we go one step further and synthesize emis- 
sion line spectra from such ab-initio coronal models (Fig. 
^ and compute average properties of these spectra, es- 
pecially emission measures and Doppler shifts (Figs. 13 
EJ, which can be considered as a test for the energetics 
and the dynamics in the model corona, respectively. This 
allows a robust comparison of the model not only to the 
Sun but also to other stars. 

2. MHD MODEL AND SPECTRAL SYNTHESIS 

The solar atmosphere is modeled from the photosphere 
to the lower corona using a sixth order fully compressible 
3D magneto-hydrodynamics (MHD) code on a staggered 
non-equidistant mesh. The simulated volume is 60 x 60 
Mm^ horizon tally and 37 Mm vertically. It includes the 
iSpitzed 1)19561) heat conductivity along the magnetic field, 
and a cooling function representing radiative losses in 
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Fig. 1. — Spatial maps in Doppler shift and intensity for the emission lines of C IV at 154. 8nm (top row) and Mg X at 62.5nm (bottom 
row). Panels (a) show Doppler shifts of the synthesized spectra as seen from straight above, panels (b) show the same for line intensity. 
This corresponds to the appearance near the center of the stellar disk. Panels (c) and (d) show side views of the computational box along 
the X and y axis in line intensity, corresponding to limb observations. The intensities / are scaled according to the average intensity (/) 
of the respective map. Please note the different color tables, especially that the C IV intensity shows much higher contrast than Mg X. 
Despite the pretty smooth appearance in intensity, the Doppler shift map of Mg X is highly structured. 



tlie optically thin corona. In the optically thick photo- 
sphere and chromosphere the temperature is kept near 
a prescribed temperature profile by a Newtonian tem- 
perature relaxation scheme. The initial magnetic field 
is a potential extrapolation of an MDI/SOHO high res- 
olution magnetogram of AR 9114 scaled to fit into the 
computational domain. The lower boundary is stressed 
by a time dependent v elocit y field, constructed from a 
Voronoi-tesselation f|Okabe et al.iil992 i ). shown to reoro- 
duce the granulation pattern ijSchriiver eral]ll997j) . The 
velocity field reproduces the geometrical pattern as well 
as the amplitude power spectrum of the velocity and the 
vorticity leaving no free parameters. 

The braiding of the magnetic field by the photo- 
spheric motions rapidly produces an intermittent corona 
in both time and space with a typical temperature of 
one million K, during the whole simulated timespan 
of ~50 minutes. The time and space averaged heat- 
ing function decreases exponentially with height, pro- 
ducing a heat input to the corona of 2— 8xlO^Wm~^ 
in agreement with coro nal energy losse s estimated from 
observations f'Withbro e & Novea 119771) . The amount 
of heating produced is constant or if anything rises 
as numerical resolution increases (Hendrix et al. 1996; 
iGalsgaard fc Nordlu nd 1996) and therefore the amount 
of heating produced in this simulation is a lower limit. 
Thus the spatial resolution of the model (400 km) is not 
sufficient to resolve the reconnection process, of course, 
but based on the above work the amount of energy input 
into the corona on larger scales should be of the right 



order. Details on and results from this MHD model can 
be found in Gudiksen & Nordhmd (2004b a). 

Using density, velocity and temperature from the MHD 
model the emissivity for a number of emission lines is 
evaluated at each grid point in the computational do- 
main. By integrating along a line of sight this finally 
results in maps of spectra when looking at the computa- 
tional box from, e.g., straight above. 

Almost all emission lines from the low corona and 
transition region are in the extreme ultraviolet (EUV; 
«500-1500A), optically thin and excited predominantly 
through electron collisions. The emissivity e is given 
through the energy hi^^ the Einstein coefficient A21 and 
the upper level population n2 for the transition from 
an upper level 2 to a lower level 1, i.e. e = /iz^ 712^21- 
To calculate the upper level density one has to evalu- 
ate the balances for excitation and ionization. In the 
latter case we assume ionization equilibrium. We have 
checked ionization and advection times and found that 
for the present models ionization equilibrium is a rea- 
sonable choice. Furthermore we assume constant abun- 
dances with photospheric values. All these quantities, 
i.e. ionization, excitation, abundances, Einstein coeffi- 
cients and transition energies, have b een evaluated us- 
ing the atomic da ta package Chianti ijDere et alJll997t 
lYoung et al.ll200^ . finally resulting in the emissivity at 
each grid point depending on the density and tempera- 
ture there. The spectral profiles at each grid point are 
assumed to be a Gaussian with a line width correspond- 
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Fig. 2. — Comparison of synthesized and observed Doppler shifts. 
The diamonds show the Doppler shifts of the temporally and spa- 
tially averaged spectra for a number of emission lines formed over 
temperatures from the transition region and low corona. The bars 
indicate the standard deviation (scatter) of the Doppler shifts in 
a time-averaged spatial map when looking from straight above at 
the computational box (cf. panels a in Fig. The thick dashed 
line shows the trend as found in observations. 
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Fig. 3. — Differential emission measure (DEM) as following from 
the synthesized line intensities compared to observations. The solid 
line shows the fit from the DEM inversion based on the spatially 
and temporally averaged synthetic spectra. The bars indicate the 
standard deviation of line intensities in a spatial map. The thick 
dashed line is based on a DEM inversion using observed quiet 
Sun disk center line radiances from SUMER IWilhelm et al.ll998l) 
scaled by a factor of two to match the active region model. 



ing to the thermal width {2k^T /m)'^/'^, where T is the 
temperature at the respective grid point and m the mass 
of the ion emitting the line. These Gaussians are shifted 
by the component of the velocity along the line of sight at 
the respective grid point. Finally we integrate the spec- 
tra along the line of sight (e.g. the vertical direction) to 
construct a map of spectra. 

These synthesized spectral maps can the be analyzed 
exactly as it is done for observations, e.g. for a spa- 
tial scan u sing the SU MER instrument on-board SOHO 
(^ilhclm et alJll995r . The resulting spatial maps in line 
intensity and line shift for a snapshot of the model run 
can be seen in panels (a) and (b) of Fig.Qlfor the lines of 
C IV at 154.8 nm and Mg X at 62.5 nm, formed around 
10^ K and 10^ K respectively. The panels (c) and (d) 
show the intensity images when viewing the box from 
the sides, i.e. along the x- and y-axis. 

3. DISCUSSION 

The synthesized corona as shown in Fig. ^ is made 
up by numerous individual loops, and the appearance is 
remarkably different in the two lines. At high tempera- 
tures the structures are rather washed out and less sharp, 
partly reflecting that heat conduction (ocT^/^) is very ef- 
ficient at high temperatures. These findings agree well 
with real observations, and the synthetic maps roughly 
resemble maps of solar observations. 

Furthermore these images reveal that in the model 
corona the low temperature emission originates not only 
from the low parts of the high temperature loops, but 
that a multitude of low lying cool dense loops are present. 
Defining the footpoints of the hot loops by circles (cov- 
ering about 2% of the area) we estimated that only some 
10 to 15% of the emission from typical transition region 
lines formed below 200 000 K originates from that area. 
Thus the transition region emission is enhanced in the 
footpoint regions of hot loops, but the majority is emit- 
ted in cooler structures. This s upports earlier sketches of 
the structure of the low corona ijDowdv et al.ll985lPeteil 



l200l . but now for the first time this multi-component 
or multi-loop structure is based on a solid model of the 
coronal structures. 

It is also important to note that the Doppler shift im- 
ages of the high-temperature line shows much more struc- 
tures than the intensity image (Fig. ^ panels a and b in 
lower row), reflecting the highly dynamic nature of the 
low corona. This shows the importance of using EUV 
spectrographs with sufficient spectral resolution to learn 
about the basic coronal processes. 

The first major test is provided through average 
Doppler shifts, which are shown in Fig. |21 The tempo- 
rally and spatially averaged shifts of synthesized spectra 
for a number of emission lines are displayed as diamonds. 
The bars represent the scatter in the spatial maps (stan- 
dard deviation; cf. Fig.^. The underlying thick dashed 
line shows the observed average variation of line shifts 
wit h temperature compi led from recent q uiet Sun stud- 
ies l|Peter fc ,Iudg3ll999D . Active regi ons ([T^ ri^^a ^IJ 
E99y) and a large number of cool stars l)Wood et alJl99'^ 
l^agano et al. 2004) show a similar behavior. The coher- 
ence of the average Doppler shifts synthesized from the 
model and the observed ones is remarkab le. Al so the 
large scatter is as found with observations jPeteil 1X999). 
For the first time the overall variation of Doppler shifts 
as a function of line formation temperature could be re- 
produced — the line shifts are caused by flows along the 
magnetic structures induced by asymetric heating. Many 
attempts have been made to understand this curve since 
the first recognition of systematic line shifts in the solar 
transition region (Doschck ct al. 197 ^. How ever, many 
of these models fafled (l|Peter fc Judggir999j) and none 
of them reached the level of qualitative and quantitative 
agreement achieved by our study. 

Only at the highest temperatures does our model 
not repro duce the recently observationally established 
blueshifts (Peter 1999: Peter fc J udge 1999). This could 
be due to the infiuence from the impenetrable upper 
boundary condition of the MHD model, which may be 
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expected to quench flows along magnetic field lines that 
intersect the upper boundary. New models, where the 
upper boundary is shifted higher up into the corona, will 
further reduce the boundary effects and could provide an 
even better match to the observed Doppler shifts. 

The other major test is to check the emission mea- 
sure distribution, which describes the emission efficiency 
at a given temperature and is defined as ir^dh/dT (for 
electron density n and temperature gradient with height 
dT/dh) fMari skalll992() . Using the synthesized lines (av- 
eraged spatially and temporally) we have performed an 
inversion to obta in the emission mea sure curve, using the 
Chianti package (lYoun g et alJ200,il) . The result is shown 
as a solid line in Fig. |21 The bars indicate the scatter of 
the emissivities this inversion is based upon (standard 
deviation). Over-plotted as a thick dashed line is an 
inversion from observed quiet Su n spectra of the sam e 
hues from SUMER fuU disk data ijWilhelm et aLllTOQ^ . 
for Mg X we evaluated disk center spectra as no full disk 
scans exist in this line. The observed quiet Sun emission 
measure is scaled by a factor of two, to be comparable 
to our active region calculation. 

The agreement of the emission measure curve from 
the observations and the synthetic spectra is remark- 
able. Previous (ID or 2D) models failed to reproduce this 
emission measure curve at low temperatures by orders of 
magnitude, basically because they had not enough emit- 
ting material below 10^ K. It has been proposed, among 
other suggestions, that in the lowest coron a a large num- 
ber of cool dense structures mig ht exist l|Gabrie]lll97^ 
iDowdv eT'al. 1986). This proposition is supported by 
our results, where enough energy is released very low in 
the corona and transition region so that enough dense 
cool structures form. Thus we can unambiguously repro- 
duce the high emission measures at low temperatures, for 
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the first time without any fine-tuning. As the emission 
measure curve is a general feature of a large variety of 
cool stars including our Sun, this suggests that the flux- 
braiding mechanism may be important in a wide variety 
of cool stars. 

Future work will have to concentrate on improving the 
coronal MHD model as well as the spectral synthesis. 
This includes the addition of the super-granular struc- 
tures to the active region patterns, which is of vital im- 
portance to understand the structure of the quiet solar 
corona. Furthermore larger scale shear motions would be 
needed to achieve a higher level of activity. In this case 
the plasma motions are expected to be faster, and there- 
fore requiring the inclusion of non-equilibrium ionization 
effects when calculating the emission line spectra. 

4. CONCLUSIONS 

This work is the first to reproduce the Doppler shifts 
and emission measure distributions in the corona of the 
Sun and many cool stars in a quantitative and quali- 
tative way. Further work will be needed to apply the 
model also to a wider variety of stellar activity, espe- 
cially when thinking of stars with different magnetic 
field structures (Donati et al. 1999) or convection pat- 
terns (l|Frevtag et al.ll2002t) on the stellar surface. This 
will show in which types of cool star coronae this heat- 
ing mechanism is the dominant one. But already we can 
conclude that flux braiding is a prime candidate for heat- 
ing in the magnetically closed parts of the coronae of the 
Sun and solar-like stars. 

Computing time for the MHD-modeling was provided 
by the Swedish National Allocations Committee and by 
the Danish Center for Scientific Computing. 
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